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We present an implementatiou of an optimal CMB trispectrum estimator which accounts for 
anisotropic noise and incomplete sky coverage. We use a general separable mode expansion which 
can and has been applied to constrain both primordial and late-time models. We validate our 
methods on large angular scales using known analytic results in the Sachs- Wolfe limit. We present 
the first near-optimal trispectrum constraints from WMAP data on the cubic term of local model 
inflation ^nl = (1.6 ± 7.0) x lif , for the equilateral model f^^'^ = (-3.11 ± 7.5) x lO" and for the 
constant model f^^^^^ = (—1.33 ± 3.62). These results, particularly the equilateral constraint, are 
relevant to a number of well-motivated models (such as DBI and K-inflation) with closely correlated 
trispectrum shapes. We also use the trispectrum signal predicted for cosmic strings to provide 
a conservative upper limit on the string tension G/x < 1.1 x 10~® (at 95% confidence), which is 
largely background and model independent. All these new trispectrum results arc consistent with 
a Gaussian Universe. We discuss the importance of constraining general classes of trispectra using 
these methods and the prospects for higher precision with the Planck satellite. 

I. INTRODUCTION 

Measurements of the cosmic microwave background (CMB) and large scale structure, such as those provided by the 
WMAP satellite or the Sloan Digital Sky Survey (SDSS), agree well with the predictions of standard single field slow- 
roll inflation. In particular the power spectrum verifies the prediction of a nearly scale invariant spectrum of adiabatic 
perturbations with a Gaussian distribution. However, there remains the prospect that significant non-Gaussianities 
may be produced by wcll-motivatived cosmological models which are consistent with the measured power spectrum. 
In order to quantify this non-Gaussianity we need to measure higher order correlators, beyond the two-point function 
or power spectrum. In [1], we presented an optimal estimation methodology to obtain the four-point correlator or 
trispectrum. This was developed from the general bispectrum estimator of [2, 3] which used separable mode expansions 
to investigate a much wider class of models than previously had been investigated, as well as to reconstruct the full 
CMB bispectrum of the Universe. These general bispectrum results were consistent with a Gaussian distribution. 
However, it is possible for the three-point correlator to remain small but for there to be a large foiir-point correlator 
(see, for example, some inflationary models [4] and cosmic strings [5]). It is our purpose here to continue to test the 
standard inflationary paradigm through an optimised and expanded search for a trispectrum signal in the WMAP 
data. 

In this paper we consider the class of four-dimensional trispectra which are independent of the diagonal. Such 
models include the cubic term of the local model, the equilateral model and the so-called 'constant' model. However, 
as we shall discuss, a general search for non-diagonal trispectrum shapes encompasses a much wider class of models, 
including most trispectra currently discussed in the literature. We demonstrate this for the non-diagonal equilateral 
shape which exhibits a high degree of correlation with the (apparently) higher dimensional shapes predicted by DBI 
and K inflation. Another example is the cosmic string trispectrum which can be reduced to a closely correlated non- 
diagonal shape. Such dimensional reduction is an important first step in testing even truly diagonal shapes because 
trispectrum estimation in the non-diagonal case is much more straightforward. Exploiting the use of a separable 
expansion, as we do here, ensures the computation is tractable, reducing the complexity from O(^niax) to 0{l'^^^^)^ 
and it also ensures the stability of the algorithms used in the analysis without the need to correct for pathological 
terms commonly present in other approaches. 

The constraints obtained here on the cubic term of the local model, the equilateral model and the 'constant' model 
result from comparison to year 5 WMAP data out to I = 500 together with a pseudo-optimal analysis of inhomogeneous 
noise and masking contributions. The estimators and other algorithms employed here were outlined in detail in [1], 
but are expressed in this paper with the simplifying assumption that the trispectrum is independent of the diagonal 
term. We validate our results by using known analytic results in the large angle limit [6] where the signal-to-noise can 
be calculated explicitly. This is important because we are able to show that previous trispectrum forecasts using this 
Sachs- Wolfe approximation were over-optimistic (see, for example, the interesting analysis of trispectrum forecasts in 
refs [7, 8]). We note that we will describe the implementation of these methods and the reconstruction of the CMB 
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trispcctrum in much more detail in a longer accompanying paper [9]. 

Some alternative approaches to extracting information about the primordial trispectrum from the WMAP data 
have been explored in the literature. In [10] the constraint —5.6 x 10^ < qnl < 6.4 x 10^ was obtained by analysing 
the A^-point probability distribution of CMB anisotropics (a non-optimally- weighted method, sec [1]). This work 
assumed a local perturbative model, $ = + fNL{'^\ + (^i)) + 9nl'^\, while using a Sachs- Wolfe approximation 
out to ^ < 1°. However, we note that our calculation of the optimal WMAP variance at Zmax ~ 500 has error bars 
of ±10.7 X 10"'' (at 95% confidence). In rcf. [11, 12], a local trispcctrum constraint -7.4 x 10^ < 5nl < 8.2 x 10^ is 
estimated from WMAP using a kurtosis power spectrum analysis (along with a constraint on tnl)- However, this 
approach does not directly subtract the effect of anisotropic noise and other systematic effects using the quadratic 
terms in the optimal trispectrum estimator [1]; we know from the present work that these are important in obtaining 
an accurate and optimized result. We will discuss and compare these alternative approaches in more detail in our 
longer paper. It is also appropriate at this stage to mention other earlier proposals for measuring the four-point 
correlator using a harmonic analysis on C0BE[13] and a wavelet approach [14]. 

n section II we review results relating primordial and CMB trispectra and their optimal estimation, here, focusing 
on primordial (and hence CMB) trispectra which are independent of the diagonal. The eigenmode decomposition of 
the trispectrum constitutes the basis of our method and is reviewed in section III. In section IV we obtain constraints 
on the cubic local model term, the equilateral model and the constant model. The accuracy of this approach is 
also verified, using large-angle analytic results for the local ^nl model. We emphasise the close correlation of our 
non-diagonal equilateral model with a broader set of equilateral-type models. We demonstrate the application of this 
method also at late times, when we constrain the cosmic string trispectrum and discuss forecasts for the Planck data. 
Finally we discuss our results and present our conclusions in section V. 



II. OPTIMAL TRISPECTRUM ESTIMATION 



CMB temperature anisotropies may be represented using the aim coefficients of a spherical harmonic decomposition 
of the cosmic microwave sky, 

-Tjri'^) = ^aimYim{i^)- (1) 

l,m 

The primordial potential $ induces the multipoles aim via a convolution with the transfer functions A;(fc), 

^AKfc)$(k)yj„(k). (2) 

The connected part of the four-point correlator of the aim gives us the trispectrum, 

Tl^m^l^m^hmsUmi = {o-l 

= {A-K)\-i)^i'' j ^^^^J — ^- Ai,{kr)/^iM^iM/^iM X 

($(ki)$(k2)$(k3)$(k4))ci;,m,(ki)y;,„,(k2)y;3„3(k3)Y;,™,(k4), (3) 

where ki = |kj| and the subscript c denotes the connected component. In this paper we will specialise to statistically 

isotropic trispectra which depend only on the wavcnumbcrs fci, fc2, fca, A:4, that is, the 'diagonal-free' trispectra class 
(and the many models which are closely correlated with these). Using this assumption we may write 

($(ki)$(k2)$(k3)$(k4))c = (27r)3 5(ki + k2 + kg + k4) Mk„k2, ^3, ^4). (4) 

Substituting this into (3), it is straightforward to show that 

J rffifi Yi,m, im^m, {n)Yism3 (A) V^m. (A) t\± , (5) 

where the 'extra'-reduced trispectrum is given by (see [6]) 

t\lu = (^) / ^^dx J {kik2k3k^f T^iktMMM) Ai^{ki)Ai^{k2)AiM^u{k4) 

X jii{kix)ji2{k2x)ji^{k3x)ji^{k4x). (6) 
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As with the bispcctrum analysis it is simpler to analyse the trispectrum - especially for scale-invariant models - in 
terms of a shape function, i.e. a scale-invariant version of the trispectrum. Removing a scaling in the diagonal-free 
case would motivate the following shape function 



S{h,k2, h, h) = {AlN)-\hk2hh)^/^ T$(fci, fc2, fcs, fc4) 



(7) 



where A<i, gives the amplitude of the Sachs- Wolfe peak of the power spectrum and iV is a normalisation factor which 
may be identified with 9tjvz,/200 such that S{k, k, k,k) = 1 (see [4]). Naively, the numerical calculation of the CMB 
trispectrum (6) appears extremely challenging (even in this diagonal- free case), as it involves a line of sight integral 
over a highly oscillatory 4D integral. However, if the shape can be represented in separable form S{ki,k2,k3,k4) = 
W{ki)X{k2)Y{k3)Z{k4), then the integral breaks down into a much simpler product of one-dimensional integrals. 

In order to compare theoretical predictions with observation it is necessary to use an estimator which sums the 
signal-to-noise over the range of multipoles probed, essentially performing a least squares fit of the data to the 
theoretical trispectrum {a]^rni^l2m2^hm3^Um4) c- [!]> the general optimal trispectrum estimator was derived giving 
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(8) 



where /sky is the sky fraction observed, the covariance matrix Ci-^mi,i2m2 = {<^iimi<ii2m2} is non-diagonal due to mode- 
mode coupling introduced by the mask and anisotropic noise, and N is the appropriate normalization N (see ref. [1] 
for further details). Issues of optimality for the bispectrum estimator corresponding to (??) were addressed in ref. [15] 
and further relevant references are reviewed in [16]. We note that a related expression for optimal estimation of the 
trispectrum in large-scale structure (and other higher-order correlators) was presented in ref. [17]. 

In this paper we follow [18, 19] by assimiing a nearly diagonal covariance matrix (Cj^mj^j,,,,-, ss 
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J and account for the noise Ni and instrument beam 6; by setting C; — > C; = bfCi + Ni 



hify^hi,^bi^t\^\^ (we drop tildes henceforth). With this identification the estimator becomes 
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where for the non-diagonal trispectrum reduces to the weighted sum 



n4 



Ci^Ci^Ci^Ci^ 



with 



{2h+l){2l2 + l)(2l3 + l){2l4 + l} 
2(47l-)2 



Pi,{lj:)PMPi,{ii)Pu{ix)dy.. (10) 



The second and third terms in (9) ensures subtraction of spurious inhomogeneous noise and masking by using the 
covariance matrix Cf™^ from an ensemble average of Gaussian maps in which these effects are incorporated. 
Again, if the theoretical trispectrum has the property of primordial separability then the summations in (9) 
and (10) become much more tractable, taking only 0{l'^^^) operations. Finally, it is important to note that the 
autocorrelator given by (10) can be used to define a natural measure on the allowed trispectrum domain and an inner 
product between bispectra {t,t') ^ Y.u ^hhhutllut'hh/'^hChCi^C^. 

It should be noted that is desirable to compare measures of the trispectra of different models. The general measure, 
tjyL, is defined in the equilateral limit for which the wavcnumbcrs fc,j and the diagonals of the quadrilateral (formed 
by the wavevectors kj) have equal values, i.e. ki = k2 = k^ = ki= ]ki -h k2] = ]ki -h ks] = k. We define 



tNL = 



9 {k f k J k'<^ k^ 



200 



(11) 



With this definition we find that for the local model t'^'jf' = 1.5tnl + ^-O^gNL (see [4]). 

As with the the shortcomings of using the parameter /jvl to compare different bispectra consistently (see [2]), the 
method of comparing trispectra using their value at a central point is problematic. Choosing a different choice of 
normalisation factor to allows us to define an integrated measure with which a consistent comparison between 
models may be made. In [6] the following choice of normalisation factor was adopted N"^ = NtNtiocB, where locB 
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refers to the g^^ expression of the local model trispectrum (with g^^ = 1 and t^vl = 0). The gpfL term is chosen 
since is independent of the diagonal term and therefore may be determined more quickly than the normalisation for 
the TjvL term. Hence we adopt the following definition for a general meaure of the trispectrum 

^NL = -^S. (12) 

This choice removes the large disparities that can exist between quoted constraints for different models. It is also 
applicable to non-scale invariant models, as well as trispectra induced hy latc-time processes such as gravitatioal 
lensing and cosmic strings. An approximation scheme for evaluating Nt using primordial shape correlators has been 
outlined in [6]. 

III. SEPARABLE MODE EXPANSION 

In [1] separable mode expansions were derived for a general (reduced) trispectrum. In this paper we consider 
trispectra which are independent of the diagonal, so the complexity in finding a general separable expansion reduces 
from O(Z^ax) to C(^max)- Here we will consider separable mode expansions at both early times and late times in order 
to project forward an arbitrary non-diagonal primordial trispectrum T(A;i, ^2, ^3, A4) and convert it into a weighted 
expansion for the late-time CMB trispectrum tj^'^. We make a separable mode expansion of the primordial shape 
function (7) in the form 

S{ki,k2,k3,k4) = Y2'^n Qn{kl,k2,k3,k4), (13) 

n 

where the Qn are product functions Qn{kx,k2,kz,k4) = <l(p{ki)qr{k2)qs{kz)cit){kA) and n = {p,r,s,t}, with (prst) 
representing the 24 cyclic permutations reflecting the underlying trispectrum symmetry. As we have emphasised for 
bispectrum mode expansions [2], the onc-dimcnsional basis functions can be chosen for convenience provided they 
are well-behaved and nearly scale-invariant. Here, we choose analogs of Legendre polynomials but with a weight 
function set by the 4D integration domain allowed by the quadrilateral condition, that is, k\ < k2 + k^ + k^ for 
ki > fc2, fca, ^4 + cyclicperms (denoted Vt)- Wc may integrate arbitrary functions /(fci, fc2, fcs, ^4) and g{ki, fc2, k^, k^) 
over this domain, defining an inner product {f,g) = /y^ f guidVT, where w is a given weight function (here to = 1). 

The Qn modes though independent and separable are not, in general, orthonormal, i.e. {Qn,Qp) = Inp 7^ ^np, so 
it is often more useful to work in terms of ortlionormal basis functions with {Rm -^p) — ^np (crccitsd from tli6 Qn 
via a Gram-Schmidt process). The relation between the Rn and Qn modes is given by 

n 

Rn = '^KpQp, where {Xnp)~^ = {RnQp) , (14) 

p=0 

where the matrices 7„p and A„p are related by 7"^^ = Xm^rp- Just as in (13), the shape function may be expanded 
in terms of the orthonormal i?„ basis S{k\,k2,kz, k^) = Ylm '^n^n{ki,k2, k^, ^4), where the = S) are related 
to the through = J2p XpnCtp. 

A similar decomposition may be applied to the late-time trispectrum f'^^ with separable product and orthonormal 

modes, denoted with a bar as Qn{h,l2,h,h) = Q{p{h)qrih)Qsih)''lt}ih) and B.nih,l2,h,h) respectively (with the 
qp{l) defined on the allowed multiple domain and a chosen ordering n O {prst}). Here, we use the late-time inner 
product defined below (10) with — We expand the trispectrum signal-to-noise in the estimator as 



<±=Y.'^Qn{hM,hh). (15) 



where, for the sake of scale invariance, we use the freedom to multiply by a separable weight function, here with 
vi = {21 + 1)"'^/*. Given the form of the (diagonal-free) primordial shape decomposition (13) we can project forward 
the Qn{ki,k2,k3,k4) modes to late times and express these in terms of the QnihjhJaU) used in (15). Explicitly, 
defining the convolved basis functions as qp{x) = f / dkqp{k) Ai{k)ji{kx) and using (6) and (13), we obtain 

c^n=^^npa^ with Tnp = j dxx^^nr{x)^-p^ (16) 
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where 7„p(a;) = {Q'^'^'"'' (x), Q^p{hMMM)) , Q[l''''''{x) = qi'{x)q','{x)q[nx)q'^{x) and q[{x) - viq[{x)/^Ci (see 
also [2]). The huge efficiency advantage of using the late-time decomposition (15) is that the line-of-sight integration 
is captured once and for all in the transformation matrix r„p and need not be repeated in any subsequent Fisher 
matrix and estimator analysis. 

Substituting the diagonal-free mode expansion (13) into (9), and exploiting (16), the estimator separates and reduces 
to the simple summation 

I^lN^_n^Q A3 AT 



E^^^n = ^H^n^^^r^n (17) 



^ n ^ n,p 

where the observed coefficients arise from separable products of the CMB map filtered with the polynomial basis 
functions qp{l), 

Explicitly, there are quartic, quadratic and constant (unconnected) parts = p''^^ — 6/3^^^ + defined by 

Pi^'> = f d^h M{p Mr Mtj , /^i^) = f cfh (A/|p A/,) Ms M,} and (3^°^ = [ dPn (M^^ M,)(M, M^}) , (19) 



where the angled brackets represent map products averaged over many Monte Carlo simulations incorporating realistic 
anisotropic noise, beam and masking effects. Remarkably, once the mode coefficients and transformations have been 
calculated, trispectrum estimation using (17) collapses down to only 0(/max) operations. 



IV. WMAP CONSTRAINTS ON NON-DIAGONAL MODELS 



We have applied the full estimator (9) using separable mode expansions (17) to the WMAP data for a variety of 

non-diagonal models. At present, this work represents the first near-optimal CMB trispectriim analysis, incorporating 
the subtraction of systematic anisotropic effects, but it neither exhausts the full precision available in the WMAP 
data nor the broad array of models which can be constrained through this general analysis. Here we report the first 
step in implementing this formalism in a CMB trispectrum analysis for examples of the simplest models (at both 
primordial and late times). We note that we employ the pseudo-optimal weighting in (9), rather than the inverse 
covariance weighting in (8), analysing the WMAP 5-year dataset out to multipoles with Zmax < 500. We coadded 
the V and W band data, using the same weights, noise model, beams and KQ75 mask as in the WMAP5 analysis. 
A similar modal bispectrum analysis was performed successfully in refs. [2, 3] and we shall describe the extensive 
work undertaken with this trispectrum estimator more fully in our longer paper [9]. However, here we note that the 

extraction of coefficients (17) from 400 Gaussian maps in a fully WMAP- realistic context produced no significant 
bias in our trispectrum estimator. Working at /max = 500, we found it sufficient to use only about fifty modes in our 
separable decompositions (13) and (15) to adequately describe all the models under investigation (i.e. rimax = 50). 



A. Local Model (cubic ^nl term) 

Local non-Gaussianity is one of the best-motivated models for non-Gaussianity because it can be produced by local 
interactions in minimal extensions of the standard model of inflation, such as those involving multiple scalar fields 
(see, for example, the comprehensive review in [20]). The gNL term of the local model arises as the coefficient of the 
cubic term in a Taylor expansion around the linear fluctuations $l: 

* = $L + ./nl($l + (*L» + ,9NL$L + (20) 

Note that there is also a local trispectrum contribution from the square of the quadratic term (in single field inflation 
this is simply tnl = (6/nl/5^), but this is much difficult to analyse because it involves integration over two solid 
angles (a treatment we shall discuss in [9]). The cubic ^nl term induces the following non-diagonal trispectrum shape 
(see [21]): 
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In the Sachs- Wolfe approximation where we take the simple transfer function A/(fc) = |ji((To — Tdoc)^), we obtain 
the large-angle solution [1] 

This analytic solution provides an excellent benchmark against which to test the accuracy of our modal decomposition 
methodology. In Figure 1 we plot the equal-/ trispectrum t\ \ for the analytic Sachs- Wolfe solution (22) and the actual 
trispectrum calculated using the full transfer functions starting from the primordial decomposition coefficients in 
(13) through the transformation matrix r„p in (15). The agreement is excellent, despite the rather pathological local 
3nl shape with divergences on faces at each ki — )■ 0. We were able to obtain an accurate primordial decomposition 
(to within a 90% correlation using a CMB-motivated cut-off as in [22]), while the r„p conversion to the late time 
decomposition (15) introduced insignificant further correlation errors (at around the 1% level). Both factors can be 
mitigated further by increasing the mode number rimax- 
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Figure 1: Plot of the ratio of the full model decomposition trispectrum with the analytic Sachs- Wolfe result (22) for the local 
qnl model. The excellent agreement for I < 60 demonstrates the accuracy of the modal formalism adopted in this paper. 



It is instructive to note the behaviour of the trispectrum in Figure 1 due to the extra power of the transfer function 
which enhances Silk damping suppression. The rapid asymptotic fall-off seen at I w 400 is actually preceded by 
a decline in the signal to noise (relative to the Ci's) which starts at / > 60 in the first acoustic peak. Making 
projections beyond / > 60 using the Sachs- Wolfe signal-to-noise estimate leads to over-optimistic expectations [7, 8], 
also indicating that the results in [12, 23] are somewhat super-optimal. The optimal (Fisher) bound at la significance 
is given by 



= 35X10'. (23) 

For our WMAP5 ^nl analysis we extracted the observed trispectrum coefficients using the estimator (17). This 

entailed subtraction of unconnected and systematic anisotropic contributions using the expectation of the filtered 

maps in (19) averaged over very many Gaussian realisations (see longer paper [9]). We note, here, that the effect 

of the anisotropic noise and mask made contributions at least at the level of the optimal variance, necessitating its 

^ 

careful subtraction (achieved here for the first time). A direct comparison of the observed orthonormal /?„ and the 
predicted orthornormal ^nl coefficients demonstrated no significant correlation between modes. Substitution into 
the modal estimator (17) yielded the constraint 

5NL = Gnl = (1.62 ± 6.98) xlO^ (24) 

The variance here was obtained from 400 Gaussian WMAP5 realisations, noting that we obtain a variance approx- 
imately 30% above the optimal bound at this resolution (in line with a previous modal bispectrum analysis [3]). 
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In order to achieve a result closer to the optimal bound the same modal methodology can be used but after pre- 
conditioning the CMB data with inverse covariance matrix weighting, as in (8). We conclude that there appears to 
be no evidence for a local gfNL non-Gaussian signal in the WMAP data. 



B. Constant Model 



The constant model represents the simplest possible primordial shape with 

fc2, ks, k4) = ^^^^^^Tconst(A;i, fc2, fcs, k4) = l. (25) 

This shape yields a CMB trispectrum t\^J^^ with features entirely due to the transfer functions. As for the bispectrum, 
it provides a useful benchmark for identifying features which can be expected to be shared by all primordial models. 
The trivial primordial decomposition with only 7^ is 100% accurate, but it still achieves above 99% accuracy for 
the late-time decomposition (15) after the Tnp transformation using 50 modes. The optimal (Fisher) bound for the 
constant model is given by 

AC£ot,t = = 1.08—^^^ = 2.67 X 10^ (26) 

Again a direct comparison of the observed WMAP spectrum /3„ and the predicted constant model indicated no 
obvious correlation. Of some interest was a negative zeroth WMAP mode /3q about which we shall comment later. 
The modal estimator (17) yielded the (equivalent) constant model constraints 

t<^^^^ = (-1.33 ± 3.62) X 10^ ^ Gnl = (-2.64 ± 7.20) x 10^ (27) 

where the t^L parameter was defined in (11) and the universally normalized trispectrum parameter Gnl was defined 
in (12). The estimated variance was again determined from 400 Gaussian simulations in the same WMAP-realistic 
context. We conclude that there is no evidence from the WMAP data for a constant primordial trispectrum signal. 



C. Equilateral Model 

Equilateral-type models are produced through the amplification of nonlinear effects around the time the modes exit 
the horizon. A non-standard kinetic term allows for such a possibility. In [4] it was shown that the leading order 
trispectrum for single field inflation models is described by a combination of three scalar-exchange trispectra, Tsi^s2,s3, 
and three contact- interaction trispectra, Tci^c2,c3 (detailed formulae for these models can be found in ref. [4]). These 
trispectra generally have a diagonal dependence, however, there is one diagonal-free model which we shall denote 
the (canonical) equilateral model with a shape function defined by 



S^{kuk2,ks,ki) = 



24A| 



TcAki,k2,k3,k4) 



{k.k^k^h)^/* 



((/ci + /C2 + fca + /c4)/4)' 



(28) 



Therefore, we shall consider the correlation of this non-diagonal equilateral model T^^ with those that have some 
diagonal dependence. This requires a full 5D Monte Carlo integration to determine the equilateral cross-correlations 
yielding (see [9] for further details): 



(Model A, Model B) 


(ci,si) 


(C1,S2) 


(C1,S3) 


(C1,C2) 


(C1,C3) 


C(Model A, Model B) 


0.95 


0.88 


0.82 


0.70 


0.56 



In particular, we note that the equilateral model correlation with the if- inflation prediction [24-26] {Tf. oc Tg^ ) is thus 
82%, while DBI inflation [4, 27] which is a linear combination of and T^^ ^2,83 exhibits an 83% correlation. Indeed, 
almost all single field inflation models of interest are well approximated by considering the equilateral T^^ model. 
Hence, we draw the important conclusion that any limit on the non-diagonal equilateral model will also strongly 
constrain other well-motivated equilateral-type models. 
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The primordial eigenmode expansion (13) with 50 modes correlates with the shape (28) at about ^ 99.7% accnracy. 
We then project this primordial shape forward to make a CMB prediction with the late-time decomposition in (16). 
A Fisher matrix analysis reveals an optimal variance for the equilateral model of 

A4topt= 5-52x106. (29) 

A direct comparison between the equilateral model and recovered mode coefficients for the WMAP5 data reveals 
little correlation, just as with the constant model. The modal estimator (17) produces the (equivalent) equilateral 
constraints 

= (-3.11 ± 7.5) X 10^ <^ Gnl = (-3.02 ± 7.27) x 10^ (30) 

We conclude that WMAP provides no evidence in favour of a wide range of equilateral-type models, including K- 
inflation, DBI inflation and all the other shapes in the table above identified above with single field inflation. 



D. Cosmic strings 

A signiflcant advantage of the formalism developed here (as in [2, 6]) is that it may be readily applied to late-time 
models as well as primordial models. Here, we apply the methodology to the trispectrum of cosmic strings. Cosmic 
strings are line-like discontinuities which may be formed during a phase transition in the very early universe [28] or 
at the end of brane inflation (in which case they are often denoted cosmic superstrings [29]). Cosmic strings are 
characterised by their tension G^. Our intention, therefore, is to use the trispectrum induced by cosmic strings to 
constrain G/U. 

In [30] the trispectrum induced by cosmic strings was derived in a WMAP and Planck context (see also [31]). The 

analysis assumed that the temperature discontinuity produced by cosmic strings is given entirely by the Gott-Kaiser- 
Stebbins effect [32, 33]. The quantity derived originally [30] was inclusive of the unconnected term and included a 
diagonal dependence. However, an accurate approximation, which correlates with the full trispectrum within 90% 
is given by 

{hhhU) v,^u -(8-G/.) ^i^i^y^yf. (0,63 ^ ^^1) 1- 2 j(,i + 500//„j ' ^^'^ 

where we note that Im = min(500,li), ^ = 1/lm and r]o/r]iss ~ 50. The values of the parameters v'^ and are found 
by simulations [34-36] to be given by the numerical values v"^ = 0.365 and = 0.42. The full trispectrum is then 
given by the permutations — pj^^^ "''^'12(3 f^^^^er details see [37]). We then (kx'omposc the cosmic 

string trispectrum in the form (15), obtaining an expansion which is found to have a 99% correlation at 50 modes. 
The signal to noise S/N = NTy/Uj2A at ;„,ax < 500 gives the following optimal bound on the cosmic string tension 
achievable using the CMB trispectrum, 

Gn < 8.8 X 10"^. (32) 

This should be compared to current CMB power spectrum constraints on Abelian-Higgs strings 7 x 10~^ [38] and on 
Nambu-Goto strings 2.5 x 10~^ [39]. We have also obtained forecasts for the optimal tension that may be probed 
using multipoles up to Zmax = 2000 (achievable using Planck data). We find for Zmax = 1000, 1500 and 2000 that the 
optimal tension probed is given respectively by G/iopt = 4.3 x 10""^, 2.8 x 10~^ and 1.8 x 10~^. Unlike the local 
gNL model the signal to noise for cosmic strings is unaffected by Silk damping. Thus, as a test for cosmic strings, we 
expect the trispectrum to provide a competitive probe to the power spectrum given the increased resolution of the 
Planck satellite and it may prove to be the most stringent probe in future surveys. 

We can compare the cosmic string trispectrum decomposition from (31) to that obtained from the WMAP5 data. 
The modal estimator (17), yields the Icr bound 

^^)' = -852±870. (33) 

where the variance is again obtained from 400 Gaussian simulations. Since the string tension must be positive, we 
deduce from this anti-correlation that cosmic strings are disfavoured by current CMB data. The 2a bound on cosmic 
strings is then given by 1.1 x 10"^. However, we note the caveat that the cosmic string spectrum is dominated by a 
positive constant term cxq > 0, which is observed in the WMAP data to be negative /3q < 0. As discussed previously 

^ 

for the bispectrum [2], the constant term /3q is susceptible to contamination by poorly resolved point sources, a 
possibility which we will continue to investigate. 
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V. DISCUSSION 

We have presented results from an implementation of an optimal CMB trispectrum estimator which employs 
separable mode expansion applicable to a wide class of isotropic inodc;ls that are diagonal independent, i.e. that 
only depend on the wavenumbers /ci, A;2, fcs, ^4. We have obtained constraints on the cubic term for the local model, 
the constant model and a notable new constraint on equilateral-type models. We found no evidence for primordial 
non-Gaussianity for these trispectrum shapes (or models closely correlated with them) at the 95% confidence level. 
The constraints on the parameter gjvzf' represent a constraint on the self-intc!raction term of the local model and, 
in addition, allows for a qualitative classification of models of the local type [40] . The constraints presented here on 
the equilateral model arc entirely new results. The importance of finding bounds on such models is, not least, as a 
consequence of the high correlation between the equilateral model, K-inflation and DBI inflation. 

In addition we have demonstrated that this estimator methodology can be applied to late-time models by searching 
for the characteristic trispectrum shape induced by cosmic strings. We have also obtained a conservative bound on 
cosmic strings G/i < 1.1 x 10~^(at 95% confidence). Using forecasts for the signal to noise at Planck resolution we have 
established that the trispectrum of cosmic strings is expected to give comparable constraints to the power spectrum 
in the near future. This constraint has the advantage of not exhibiting degeneracies with background cosmological 
parameters, while also being largely independent of the underlying properties of the cosmic strings (unlike gravitational 
wave constraints). 

This paper marks a significant first step in a general analysis of trispectrum models. Among the advantages of the 
modal approach pursued here are the exploitation of efficiencies arising from separability and the transfer function 
transformation, the absence of pathologies when representing shapes, and the opportunity to investigate non-separable 
models which were previously deemed intractable. We also note that this general mode expansion allows for an 
accurate characterisation of the noise and foregrounds which must be subtracted to achieve an optimal estimation 
measure. While the constraints published in this paper are consistent with Gaussianity, the extension of the analysis to 
general trispectra represents an important development [9]. Such an analysis will, for instance, allow a measurement 
of the local trispectrum parameter tjvl providing a test of local inflation which requires tjvl > (6/jvl/5)^. An 
implementation of the late-time modal estimator outlined in [6] to include trispectra dependent on the diagonal term 
will allow identification of any trispectrum whether generated at primordial times like inflation or late-times like 
gravitational lensing or second-order gravitational effects. This, in conjunction with recent results classifying CMB 
bispectrum constraints [3], offers the hope of a comprehensive test for non-Gaussianity. 

We conclude by noting that we have obtained near-optimal trispectrum constraints on primordial local gNL and 
equilateral (and correlated) models, together with late-time constraints on cosmic strings. We find no evidence for a 
significant trispectrum in these non-diagonal cases, about which we will provide greater detail elsewhere [9], together 
with a reconstruction of the CMB trispectrum obtained from the WMAP data. The new results here represent a 
significant test of the simplest standard model of inflation, further affirming the Gaussian hypothesis for primordial 
fluctuations. 
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